Bulk and Spectral Observables in Lattice QCD 
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Abstract. We review recent developments in lattice siumulations of the equation 
of state, order of the thermal phase transition and the determination of the pseudo- 
critical temperature in (2+l)-flavor QCD. Owing to the increasing computer power, 
new argothithms, and improved fermion formulations, studies of bulk QCD matter are 
approaching to the stage of precision science. We also review recent lattice studies on 
the spectral properties of heavy quarkoniums inside the quark-gluon plasma (QGP). 
Although they are still in an exploratory stage, interesting physics in relation to the 
strongly correlated QGP is being extracted. 



1. Introduction 

One of the main goals of the lattice QCD studies is to make first principle analysis 
of hot/dense QCD and to supply reliable inputs to quark-gluon plasma (QGP) 
phenomenologies pp. Shown in Fig{T] is a schematic QCD phase diagram indicating 
three basic phases; the chiral symmetry broken (xSB) phase, the color superconducting 
(CSC) phase and the QGP phase. Precise locations of the phase boundaries and the 
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Figure 1. A schematic QCD phase diagram in the plane of temperature T and baryon 
chemical potential ^t B . Solid lines indicate the first order phase boundaries. Several 
critical points at which the first order lines terminate may exit, e.g. the high T critical 
point [2] and high fj, B critical point [3] . 
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critical points as well as the dynamics in each phase should be determined by non- 
perturbative method such as the lattice QCD. 

QCD partition function in a finite box with a spatial volume V and the lattice 
spacing a can be written as 



where U is the gauge field defined as an element of SU(3) C and F(U) denotes the quark 
contribution to Z. Eventually, we need to take the continuum limit (a — > 0) and the 
thermodynamic limit (V — > oo) to extract physical observables in the real world. 

In the past few years, there have been considerable progress in lattice QCD 
approach: (i) available computer speed becomes as fast as 50 Tflops, (ii) improved 
fermion formulations have been tested extensively (such as staut, asqtad and p4 
improved actions for staggered quarks and the clover action for Wilson quark), (iii) 
new fermion formulations with good chiral properties are started to be used (such as 
domain wall and overlap fermions), and (iv) new algorithms for full QCD simulations are 
proposed and implemented (such as the rational hybrid Monte Carlo method (RHMC) 
and the domain-decomposition hybrid Monte Carlo method (DDHMC) |5J). Shown in 
Figf5]is so called the Berlin wall plot where the computational cost in unit of Tflops-year 
is plotted against m^/nip. The quark mass dependence of the cost becomes weaker than 
before, and simulations with realistic quark masses may be done in a reasonable amount 
of time by using O(50) Tflops machines. 

In the following, we will focus on the bulk and spectral properties associated with 
the QCD phase transition at finite T with /x B = 0. For the developments in lattice QCD 
at finite chemical potential, see the recent review 0j. 
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Figure 2. Updated Berlin wall plot which shows the cost to generate 1000 independent 
gauge configurations on 24 3 x 40 lattice with a = 0.08 fm in various fermion 
formulations with modern algorithms [5] . 
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The equation of state at finite T characterizes the bulk QCD matter. Also it becomes 
an fundamental input together with the transport coefficients [7] to the relativistic 
hydrodynamics for QCD fluid \8\. Shown in FigJH] is the state-fo-the-art calculation 
of the energy density e as a function of T in (2+1) -flavor QCD with asktad improved 
staggered quarks. (See also the recent results for 3-flavor QCD with the p4 action [10] 
and for (2 + l)-flavor QCD with the stout action |llj.) 

The figure shows that ejT A increases rapidly at T ~ 200 MeV and approaches the 
Stefan-Boltzmann (SB) limit from below. The deviation from the SB limit at high T is 
due to the quark-gluon interaction which decreases only logarithmically as T increases. 
Quantitative understanding of the deviation is a non-trivial issue and cannot be treated 
in naive thermal perturbation [13J. 



3. Order of the transition in (2+1) -flavor QCD 

Using the pressure P(K) = (T/V) ■ \nZ and a vector K = (T, /i, m q , a, V, ■ ■ •), one 
may define the n-th order transition as the one where non-analyticity of P in the 
thermodynamic limit (V — > oo) first appears in the n-th derivative of P with respect 
to K. For example, if (d/dT)P = s is discontinuous at certain T, it is the first order 
transition. On the other hand, if (d/dT)P is continuous and (d/dT) 2 P = (d/dT)s = 
c v /T diverges at certain T, it is an example of the second order transition. If smooth 
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Figure 3. Energy density e as a function of T for (2+l)-flavor QCD with asqtad 
improved action on the lattice, N t x = 4 x (12 3 - 16 3 ) and 6 x (12 3 - 20 3 ) [9]. 
The strange quark mass m s is about the physical value, while the light quark mass 
m u d — 0.1m s corresponding to m n /m p ~ 0.3. The Sommer scale n = 0.138(7)(4) 
fm is used to set the scale. Lattice spacing at T ~ 200 MeV is 0.247 (0.165) fm for 
N t = 4 (6). 
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phase change takes place without the non-analyticity of P, it is the crossover. To 
make a distinction between the phase transition and the crossover from the numerical 
simulations, the scaling behavior of certain observables as a function of V needs to be 
studied (the finite size scaling analysis) [IT]. In the (2+l)-flavor QCD, the relevant 
observable is the chiral susceptibility whose volume dependence at its maximum reads 

Xm = (d / dm uA f P -> V a at the peak, (2) 

where a ~ 1 (2/3) for the first (second) order transition, and a ~ for crossover. 

Shown in the left panel of Fig H] is the dimensionless chiral susceptibility Xm/T 2 
around its peak for different values of the dimensionless volume N s = V/a 3 = 
18 3 , 24 3 , 32 3 for (2+l)-flavor QCD with realistic masses of u, d, s quarks |12j . There is 
no appreciable volume dependence for the height of the susceptibility, which suggests 
the crossover nature. The right panel of Figj4] shows the subtracted chiral susceptibility 
in the continuum limit Ax ( = the peak of % m (T) — x m (0)) in a renormalization group 
invariant combination T 4 / (m^ d Ax) as a function of the dimensionless inverse- volume, 
1/(T 3 V). The height of the susceptibility approaches to a non- vanishing constant in 
the thermodynamic limit and has different scaling behavior from those of first or second 
order transitions. This gives a strong evidence that the hadron-QGP transition at finite 
T with zero chemical potential is a crossover in the real world. 
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Figure 4. Left: The chiral susceptibility in (2+l)-flavor QCD as a function of the bare 
coupling, 6/g 2 , for different values of the spatial volume, 7V t 3 = V/a 3 = 18 3 ,24 3 ,32 3 . 
The stout improved staggered fermion action is employed and the quark masses 
are chosen to reproduce m K /m n ~3.689 and / K /m w =1.185. Right: The volume 
dependence of the peak of the subtracted chiral susceptibility in the continuum limit. 
Dotted (dashed) line is an expectation from the first order (second order) phase 
transition. Figures are taken from [12] . 
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4. Pseudo-critical temperature for (2+l)-flavor QCD 

Since the hadron-QGP transition in the real world is crossover, the "critical" 
temperature T c is not a well-defined concept. Nevertheless, one may introduce a 
"pseudo-critical" temperature peak position of certain susceptibility. The actual 

value of T pc is different for different choice of susceptibility, e.g. Xm = (d/dm n< i) 2 P, 
Xt = (d/dT) 2 P, and XmT = {d 2 /dmdT)P. Note also that multiplying an arbitrary 
function of m u d and T to the susceptibilities can change the value of T pc . 

Shown in Figj5] is a summary of (pseudo-) critical temperature recently 
calculated from the dimensionless chiral susceptibility, Xm/T 2 , for (2+l)-flavor QCD. 
Extrapolation to the continuum limit is taken by using the data for several values of 
N t . The upper two bars are T c with different extrapolations to the chiral limit m U( j = 
by MILC collaboration [Hj. The middle bar is T pc extrapolated to the realistic m u d 
by RBC-Bielefeld collaboration [15]. The lower bar is T pc calculated at the realistic 
m u d by Wuppertal-Budapest collaboration [16]. Note that different improved actions 
and different scale determinations are adopted in three groups. When extracting T pc in 
MeV unit, we need to use the chain rule: 



Here M, which has a dimension of mass, is chosen either to be a direct physical 
observable (such as m p , f n and f K ) or to be a semi-empirical quantity (such as 
the Sommer scales 1/rn and 1/ri defined from the heavy-quark potential through 
r 2 dV(r)/dr\ r=ro = 1.65 and r 2 dV(r)/dr\ r = 1.0). 
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Figure 5. Recent determinations of the (pseudo-)critical temperature for (2+1)- 
flavor QCD with staggered quarks. As for the lowerst bar, we add 9 MeV to make a 
comparison to other results obtained from Xm/T 2 [IB] . 
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There are a number of issues to be studied further before drawing definite conclusion 
on the value of T pc in the real world. In particular, (i) careful continuum extrapolation 
(a — > 0) by the data with larger N t , (ii) determination of T pc in the case of Wilson fermion 
[18] , and (iii) resolving the difference between the Sommer scale in the continuum limit 
obtained from a being determined by the 2S-1S mass splitting of the bottomonium with 
staggered quark (r = 0.468(7) fm) [19] and that obtained from a being determined by 
the kaon mass with Wilson quark (r = 0.516(21) fm) [20] . 



5. Heavy quarks inside QGP 

Heavy quarks in the real world (such as the charm and bottom) have finite masses 
and may receive substantial kick from the light plasma-constituents which have typical 
energy of about 3T. Therefore, it is important to examine the dynamical correlations 
(both spatial and temporal) of the heavy quarks inside QGP. It is also relevant to study 
the modification of the heavy quarkonium properties in hot QCD matter |21j . 



5.1. Charmonium wave function in quenched QCD 

An attempt to study the J/^ wave function at finite temperature in quenched QCD 
(F(U) is take to be unity in Eq.([T])) on an anisotropic lattice was performed by 
Umeda et al. [22J. They have studied the equal-time Bethe-Salpeter wave function 
w(r, t) at finite T in the Coulomb gauge. In Figj6j the normalized wave function, 
4>(f, t) = Wi(f, r)/wi(0, r) is plotted below and above T c for various time slices. The 
black dashed line is the wave function prepared at r = 0. For free quarks without gauge 




Figure 6. The equal-time Bcthc-Salpctcr wave function normalized at the origin as a 
function of the spatial separation r between the charm and the anti-charm. Blue (red) 
points correspond to free (interacting) quarks [22] . The black dashed line is the wave 
function initially prepared at t = 0. Note that t = r in this figure. 
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interactions, the wave function becomes broader as the imaginary time r increases. On 
the other hand, the wave function for the interacting system is almost unchanged even 
at T = 1.52T C , which suggests that J/^ may survive as a bound state in the deconfined 
phase. 



5.2. Charmonium spectral function in quenched QCD 



The spectral functions of hadronic correlators give us another information on hadronic 
modes at finite T as originally suggested in [23]. The spectral function o~(u,p) for 
hadronic correlation may be defined through the spectral decomposition: 



D(r,p) 



a(u,p) duo (0 < r < T 



(4) 



The maximum entropy method (MEM) provides an efficient and powerful way to obtain 
a unique a from the lattice data D without making a priori parameterization [23] . 

Applications of MEM to the ss mesons at finite T [25] and to the charmoniums 
at finite T [26l [27] [28] [29] within the quenched QCD have been carried out and it was 
shown that the J/ty and r\ c survive even above T c as distinct peaks. Shown in the left 
(right) panel of Fig 17] is the dimensionless spectral function for J/ty (i] c ) below and above 
T c on 32 3 x (96 — 32) (24 3 x (160 — 34)) anisotropic lattices with the spatial volume 
V = (1.25fm) 3 (V = (1.34fm) 3 ) taken from Ref.[26] (Ref.pj). The figures indicate that 
the low-lying s-wave resonances still survive even at 1.5 T c , which is consistent with the 
result in FigEJ Note that this result is not an artifact of the small spatial volume, which 
was carefully checked in [30] . 

There are several topics which are interesting to be studied within the quenched 
QCD, such as the spectral functions of the charmoniums moving inside the hot medium 
[3T] and the spectral functions of bottomoniums above T c |32j . 
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Figure 7. Left: Spectral functions measured in quenched QCD simulations on an 
anisotropic lattice with the spatial volume (1.25 fm) 3 for the J r /\E f channel Right: 
Spectral functions measured in quenched QCD simulations on an anisotropic lattice 
with the spatial volume (1.34 fm) 3 for the r/ c channel [2"3] , 
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Figure 8. Spectral functions in the vector and pseudo-scalar channels in 2-fiavor QCD 
on 8 3 x (32 — 16) anisotropic lattices. The spatial lattice volume is V ~ (1.2 fm) 3 and 
the light quark mass corresponds to m^/nip ~ 0.5. The figures are adapted from 34J. 



5.3. Charmonium spectral function in 2- flavor QCD 

In quenched QCD, only gluons are thermally excited in the plasma, while, in full QCD, 
both quarks and gluons are thermally active. Because the number of degrees of freedom 
increase in full QCD, the dissociation rate of the chamoniums will be larger. On the other 
hand, the pseudo-critical temperature T pc in full QCD (160-200 MeV) is substantially 
smaller than T c ~ 270 MeV in quenched QCD. Due to the compensation of these effects, 
the ratio of the dissociation temperature T* and the (pseudo-) critical temperature in 
full QCD is not much different from the ratio in quenched QCD and could take a value 
of around 2 [33J. 

Aarts et al. has recently reported an exploratory study of 2-fiavor QCD on 
8 3 x (32 — 16) anisotropic lattices |34j . Shown in FigJH] are the spectral functions of 
the J/\I/ and rj c . The results suggest that, even with the thermal excitations of light 
quarks, the low-lying charmoniums may survive inside the quark-gluon plasma above 
T c until they dissociate away around 2T pc . 

6. Summary 

We reviewed recent developments in lattice simulations on the equation of state, 
crossover nature of the thermal phase transition and the determination of the pseudo- 
critical temperature in (2+1) -flavor QCD. Recent lattice results on the spectral 
properties of the heavy quarkoniums inside the quark-gluon plasma are also summarized. 
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